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We study a superfluid Bose system with single-particle and pair condensates on the basis of a half-phenome- 
nological theory of a Bose liquid not involving the weakness of interparticle interaction. The coupled equations 
describing the equilibrium state of such system are derived from the variational principle for entropy. These 
equations are analyzed at zero temperature both analytically and numerically. It is shown that the fraction of 
particles in the single-particle and pair condensates essentially depends on the total density of the system. At 
densities attainable in condensates of alkali-metal atoms, almost all particles are in the single-particle conden- 
sate. The pair condensate fraction grows with an increasing total density and becomes dominant. It is shown 
that at density of liquid helium, the single-particle condensate fraction is less than 10%, which agrees with ex- 
perimental data on inelastic neutron scattering, Monte Carlo calculations and other theoretical predictions. The 
ground state energy, pressure, and compressibility are found for the system under consideration. The spectrum 
of single-particle excitations is also analyzed. 
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1. Introduction 

The first experimental observations of Bose-Einstein condensation in dilute gases of allcali-metal 
atoms have stimulated a great interest to this remarkable phenomenon manifested also in super- 
fluids and superconductors. However, in spite of a significant progress |4- 8] in the study of Bose systems 
with condensate, their theory is far from being completed. Among the open theoretical problems we men- 
tion the following ones: description of Bose systems with strong interaction, microscopic justification of 
observable excitation spectrum in a superfluid ^He, role of single-particle and pair condensates in the 
phenomenon of superfluidity, the existence of elementary excitations with activation energy along with 
sound excitations. 

As was first showed by Bogolyubov |9], the interparticle interaction essentially affects the behavior 
of a many-body Bose system at low temperature. In particular, in consequence of a weak interaction, 
the number of particles in a condensate at zero temperature is not equal to the total particle number 
in the system. This effect, usually referred to as depletion of a condensate, is caused by the presence of 
pair anomalous averages which are similar to Cooper correlations in superconductors. The pair corre- 
lated bosonic atoms form a pair condensate which, along with a single-particle condensate, specifies the 
superfluid density. The role of pair correlations in superfluid Bose systems has been studied by many au- 
thors ilol - Eill . Superfluidity has also been treated in terms of pair condensation only, in the total absence 
of a single-particle condensate |20- 22]. 

Note that accounting for pair correlations in a Bose system can produce a gap in the spectrum of el- 
ementary excitations Ill [ 4l 41l. T he existence of a gap essentially depends on the used approximation or 



truncated Hamiltonian (5|,|8|,|23|,|2^. In particular, the well-known Hartree-Fock-Bogolyubov approxima- 



tion extended to bosons generates a gap in the spectrum. In virtue of the Hugenholtz-Pines theorem (25 
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according to which the spectrum should be gapless, this fact is usually considered as a defect of the men- 
tioned approximation. Therefore, some efforts have been done to reformulate the Hartree-Fock-Bogolyu- 
bov theory in order to remove a gap in the single-particle excitation spectrum (see (2^ and references 
therein). However, such schemes involve additional assumptions or parameters which does not alow 
one to treat them as a rigorous solution of the problem. Note that Pines |27] has pointed out that the 
correctness of their theorem depends on the validity of expansions into a series of perturbation theory. 
Moreover, the Hugenholtz-Pines theorem cannot be applied to all models because it is valid for a specific 
truncated Hamiltonian used by the authors. The subsequent analysis of Bogolyubov's -theorem has 
showed that no general conclusion can be made concerning the excitation spectrum of a superfluid and 
long-wavelength density excitations are insensitive to (7(1) symmetry breaking flS]. As was also stressed 
by Bogoljrubov and Bogolyubov (jr.), the presence of correlated pairs in a Bose gas of interacting parti- 
cles leads to the fact that the spectrum consists of two branches — the branch with activation energy (or 
with a gap) and a phonon mode fl^. The similar form of spectrum for a superfluid *He was also pro- 
posed within a phenomenological model |30]. Perhaps, one may consider that the existence of a gap in 
the spectrum of single-particle excitations does not contradict the general physical principles. Moreover, 
the possible existence of elementary excitations with activation energy along with sound excitations re- 
quires further theoretical and experimental investigations. The existence of a gap in the sing le-particle 
excitation spectrum has been recently discussed in terms of variational wave functions 13111 and with- 
out the involvement of a conjecture of c-number representation of creation and annihilation operators 
with zero momentum |32]. Note that the c-number representation results in a reduction of Fock space, 
though one can prove that it correctly reproduces the pressure (i^ and condensate density js^. tUl in 
thermodynamic limit. In a recent study \3ff\, it has been shown that the restoration of the state with zero 
momentum to the Hilbert space with subsequent exact numerical diagonalization of the total single-mode 
Hamiltonian yields the excitation spectrum with a finite gap. 

The description of many-body systems that does not require the weakness of interparticle interaction 
can be formulated using the quasiparticle concept. The well-known example of such a description is 
the Fermi liquid theory proposed by Landau (37|] and Silin | 38]. This theory has been extended to the 
study of various superfluid states of Fermi systems by introducing pair anomalous averages liiljilll . 
The main advantage of the extended half-phenomenological approach is that it is valid for an arbitrary 
energy functional and does not restrict the intensity of interaction. It can be applied to a wide range 
of systems, including strongly interacting superfluid nuclear matter (izl - li^l . On microscopic level, such 
half-phenomenological description is equivalent to a self-consistent mean field theory |45]. Subsequently, 
the ideas of the extended Fermi liquid approach have also been disseminated to a Bose superfluid i4lll46il 
with a corresponding microscopic justification |14]. Note that the self-consistent mean field model is an 
effective zeroth-order approximation of quantum-field perturbative theory |47]. 

In the present paper, we study a superfluid Bose system with single-particle and pair condensates 
without any restriction on intensity of interaction and conjecture of c-number representation of creation 
and annihilation operators. To this end, we employ the formalism |41, 46] developed by us earlier. From 
the principle of maximum entropy, we present a brief derivation of the coupled equations describing 
the system at finite temperature. The obtained equations are analyzed in detail both analytically and 
numerically in case of zero temperature and contact interaction. It is shown that for systems of low 
density (dilute gases of alkali-metal atoms), the pair condensate fraction is negligibly small in comparison 
with the fraction of single-particle condensate. Such systems, as expected, are well described by the Gross- 
Pitaevskii equation |4, 6, 7]. The role of pair condensate becomes determinative for dense systems. In 
particular, it is found that the single-particle condensate fraction in a superfluid *He is less than 10% 
at zero temperature. This result is in a good agreement with experiments lisUsi!) on inelastic neutron 
scattering, with Monte Carlo calculations 1 52] and with theoretical approach that involves the calculation 
of single-particle density matrix expressed through the structure factor l53-t55il. We also calculate the 
ground state energy, pressure, speed of sound for a system with two condensates. The single-particle 
excitation spectrum is analyzed. 
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2. Formalism and basic equations 

Consider the basic principles that underlie the theory of superfluid Bose systems with single-particle 
and pair condensates (for details see the reviews ^il'.Ee']). This theory is formulated in close analogy with 
extended Fermi-liquid approach to superfluids 1 39-41]. The state of a superfluid Bose system is specified 
by the condensate amplitudes, 

bp-Sppap, bp-Sppa^, (2.1) 
as well as by normal and anomalous single-particle density matrices, 

/pp' = Sppfl^,flp, gpp. = Spp<2p-ap, gpp, = Sppa^,a],, (2.2) 

where creation and annihilation operators Up, ap satisfy the usual Bose commutation relations and p is 
a statistical operator of the system which we define below by equation <2.3> . Note that fppi - f*,^ and 
gppi - gp/p. The condensate amplitudes 12. 1> describe a single-particle condensate, while anomalous av- 
erages l l2.2t characterize the pair correlations between particles and indicate the existence of a pair con- 
densate. 

We will approximate the statistical operator p by operator that contains a general quadratic form of 
flp, as well as linear terms in creation and annihilation operators, 

p = exp(Z-F), F^a^Aa+^\^aBa+a^B*a^^ + a^C+C*a. (2.3) 

Here, we have omitted the repeated summation indices bearing in mind that, e.g., a^Aa = flpApp'^p' or 
fl+C = flpCp and Z is found from the normalization condition, Spp = 1. The quantities ^pp', Bppt, Cp are 
related to /pp/, gpp^, bp by equations 12. 2> . 12. It . The terms linear in creation and annihilation operators 
that stand in the exponent ofp can be removed by the unitary transformation of a c-number shift (for 
details see references I 4ll 14^.15^ ): 

Uap [/^ = flp -H fop , Ual = flj, -H fop . (2.4) 
Then, the statistical operator g - UpU^ wiU include only quadratic terms in a^, Up as well as matrices 



g - exp 



Z - fl^ Afl — ^aBa+ a^B* a^] 



Z^Z+b*Ab+^[bBb+b*B*b*). (2.5) 



2 

For this statistical operator, we have Sppflp = 0. Therefore, according to equations 12.4> . one obtains 

Sp pflp = Sp p(flp -I- fop) = fop 7i 0. 

Using the unitary transformation 12.4) and equations 12.2) . it is easy to introduce the correlation functions 

, and e'^ ,: 
PP pp 

/pp' = Sp paj,, flp = fop, fop + /p'p, , 

gpp' = Sp pflp' flp = fop' fop + gpp, , (2.6) 

where 

/pp' = Sp pfl^, flp , gpp, = Sp pflp' flp . (2.7) 

From equations 12.7) . which relate the matrices Appi,Bppi to /pp,,gpp„ it follows (s^ thatp = p(/'^,g'^,g'^^)- 
Therefore, the entropy of a Bose system S = - Sp p In p = - Sp p In p is a functional of correlation functions 
only, S - S(/'^,g'^,g'^^). Moreover, it can be shown |41, 46] that the matrices Appi and _Bpp' represent the 
derivatives of entropy with respect to correlation functions: 

dS dS I dS I f 



Qfc PP ' g c 2 PP ' ae'^ 2 PP' 

•'p'p Sp'p "gp'p 
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Taking into account equations <2.8) . i2J\ we can see that these ntatrices are expressed through the av- 
erages of creation and annihilation operators. Therefore, they are also transformed under the global 
phase transforntations of ap, a^, so that statistical operator <2.5> is invariant. In addition, such choice 
of p (or g) is based on the fact that it satisfies the principle of spatial correlation weakening and Wick's 
theorem applies to it |56]. The introduced statistical operator defined by equation <2.5) describes a su- 
perfluid many-body system of interacting particles in the language of free particles (or quasiparticles) 
with a modified dispersion law. It is worth stressing that the efficiency of such description is essentially 
determined by the choice of coefficients ^pp^ and Bppi. As we will see below, these coefficients are chosen 
to satisfy the requirement of maximum entropy for fixed values of additive integrals of motion (or con- 
served quantities). In the language of a self-consistent mean field theory, this requirement is equivalent 
to the fact that the Hamiltonian of such approximation is the closest to the exact one |14]. Therefore, the 
statistical operator <2.5) and the maximum entropy principle give the most accurate description of the 
system within the quasiparticle approximation even if the interparticle interaction is not weak. 

A compact formulation of a theory under consideration is given in terms of a two-row matrix f'^ that 
combines the correlation functions and a vector if/, whose components are the condensate amplitudes, 

^^^[e^' -I'-r)' ^^(^')' 

where tilde denotes the transposed matrix. Thermodynamic equilibrium of the system is determined by 
the maximum of entropy S = -Sp p In p for fixed values of conserved quantities, such as energy, momen- 
tum, and total particle number. Using the unitary transformation I l2.4t and u - v transformations, one 
can show |41, 46] that the statistical operator I l2.5t is reduced to a diagonal form po- Then, for the entropy 
of the system S - -Sp polnpo, the following combinatorial expression is vaUd: 

S = -tr[/oln/o-(l + /o)ln(l + /o)], /o^Sppofl^fl, 

which allows it to be expressed through the introduced two-row matrix fillfi^l : 

S(/'^) = -ReTr/'^ln/^ (2.10) 

Here, tr... is a trace over the momentum variables which specify a single-particle state, while Tr... is 
taken over the the two-row matrix as well as the momentum variables. The introduced entropy depends 
on correlation functions only and does not depend on the condensate amplitudes. 

The energy of a superfluid Bose system is a functional of correlation functions and condensate ampli- 
tudes, E - Eif'^ yif/). It is obtained by averaging a microscopic Hamiltonian, 

E{f',ii/)^SppH[cil,ap]. (2.11) 

In fact, the normal-ordered Hamiltonian H^a^, Qpj of a Bose system can take into account binary, triple, 
and higher-order interactions of particles. However, we will use the following second quantized Hamil- 
tonian with binary interparticle interaction, 

Hial,ap\^Y.^^l'^P + TU E v^Pi-P3)«pi42'^P3'^P4^Pi+P2,P3+P4- (2.12) 
' p '^in ^ V p[...p4 

Now, we find the explicit form of the energy functional E{f'^,ifr) for the written Hamiltonian. Due to 
unitary transformation <2.4) . the averaging given by equation 1 12. 11> is reduced to 



E 4f) = SppH [al + b*,ap + fop) . 



Since g has a Gaussian form that involves both normal and anomalous pairs of second quantized opera- 
tors in the exponent, we can apply Wick's theorem with non-vanishing pairwise normal and anomalous 
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averages [or contractions, see equations (2.7H . Thus, we have t57il 

piP2 



+ ■^7 L /piP2[^(Pl-P25 + ^(P2-P4)]^p3V^P3+P2,Pl+P4 



^Pl-P4 



+ ^ E /pip2/p'3P4t'*'(Pl-P2) + '^'fP2-P3)]5pi+p3,P2 + P4 
^1' P1...P4 



+ ^ E V(P2-P3)[gpip2^p3^'p4+h.C.|5p,+p,,p3+p, 
P1-P4 

+ 777 E gpip2gp3P4['^(Pi-P3) + '^'(Pi-P4)l<5pi+P2,P3+P4. (2.13) 



where the first term i?(i/>) = Eibp, bp) is constructed from the condensate amplitudes bp, b* only, 

p2 Y 

E'^V)^lL:r-b;bp + — E v(pi -p3)fo;jfo;,fop3fop45pi+p2,p3+p4 . (2.14) 

p Z/7I Z 1/ pj...p4 

It can be proved that E(i(r] is related to the total energy functional by the differential opera- 

tor t4l„46,1. 

E{f',{ff)^RE{^r), 

where 

I d d 1 d , d Id „ d 
^[db-' db* 2db^ db 2db*^ db* 

In the expression for R, we have omitted the repeated summation indices, like in equations I l2.3t . ( 12. 5t . 

The total particle number, similar to entropy )2.10> . can be expressed through the two-row matrix f'^. 
Like in quantum mechanics, in which the physical quantities correspond to operators, we assume that a 
physical quantity a is associated with a two-row matrix a. Then, the average of a physical quantity a is 
given bv L4L.4&1 

(a) = tr/a = - (Tr/'^a- tra-i- i/'*f3ai/>) , a=| ^ ^_ 

where is the Pauli matrix and tr fa = Zp(/'2)pp- Since is the generator of unitary gauge transfor- 
mations, it should be interpreted as the particle number operator. Therefore, according to the above 
formula, we have 

W(/SV.) = tr/=i(Trff3-trl + rV'), f 3 = ( J \ 
where 1 = Sppi . The calculation of the traces gives the following expression for the total particle number: 

N{f',v)-L(fp+K^p)' /PP^/P- (2.15) 
P 

As we have already mentioned, equations that determine the equilibrium values of correlation func- 
tions and condensate amplitudes are obtained from the principle of maximum entropy for fixed values of 
the additive integrals of motion — energy, total particle number, and total momentum. However, below 
we will study the system at rest (we do not introduce the latter integral of motion). Then, the problem of 
conditional maximization of the entropy can be reduced to the problem of unconditional minimization 
of the following non-equilibrium thermodynamic potential: 

n(f,i^) = -S(^) + ;6[£(f ,t^)-pAr(/^i^)], (2.16) 
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where /3, fifi are the corresponding Lagrange multipliers (fi - l/T is the reciprocal temperature and ^ 
is the chemical potential). The solution of the formulated variational problem gives the self-consistency 
equations jililigil : 

f ^[expP{£-fi)-l]-\ (2.17) 
fj-j^iff^O, (2.18) 



where 
and 



(2.19) 



dE{r,v] jE[r,v] dE[r,<if] 

App, = 2^^, r,,^^^. (2.20) 

Equation 12.17> has a natural form — namely, it reflects the fact that the matrix f'^, which is constructed 
from the correlation functions, has a structure of the Bose distribution function if £ is interpreted as the 
operator of quasiparticle energy. Note that e and /"^ are Hermitian matrices in indefinite metrics intro- 
duced in references Jill Equation <2.18) has a structure similar to the stationary Gross-Pitaevskii 
equation without external potential. However, the principal difference is that the energy of the system 
depends now not only on the condensate amplitudes but also on correlation functions. The coupled equa- 
tions 1 12. 17L 1 12. 18> describe an inhomogeneous equilibrium state of a superfluid Bose system with single- 
particle and pair condensates if the energy functional E (/'^, ■^/] is known. 

In conclusion of this section we would like to note that a similar approach has been developed to 
extend the theory for a normal Fermi liquid to superfluid states I's^-'il]. In this case {Up) - (Up) - and 
the system is described by one equation only, which has a structure of the Fermi distribution function 
[similar to equation 12.17H . 



3. Spatially homogeneous state 

In the case of a homogeneous system, the self-consistency equations 1 12.1 7> -l [2^20) have a more simple 
form. In particular, the correlation functions and condensate amphtudes meet the following relations: 

fpp'^fp^P.p'' g'pp'^gpVp'' bp^bo5p,o, (3.1) 

where fp - fp* and gp - g£p. These relations show that only particles with zero momentum and pairs of 
particles with total zero momentum are in a condensate. Then, in accordance with equations d2.20t . we 
have 

£pp' = £p^p,p' . App; - Ap5p _p; , (3.2) 

where 

Let us transform equation <2.17) taking into account the above conditions for spatial homogeneity. To 
this end, we expand the 2x2 matrix P^-P[e- p.), entering this equation, in terms of the traceless PauU 
matrices ft and the identity matrix f q: 

^^^^Cifi + cofo, (3.4) 

where ^-E-fi. We have employed the fact that for a spatially homogeneous state the first relation from 
13.2\ yields f = ^. Taking the trace of both sides of equation I |3.4> we have cq-O and, consequently, 

[ciTif^p^Ei, Ei^ep-Ap^;- 
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Let «(x) - [e^ - 1) ^ be the Bose distribution function. Then, equation <2.17t can be written in the form 
f'^ - n ()3f). Bearing in mind equation <3.4> . we can find the following representation for this equation: 

f = i [n[(iE^] + n[-pE^]\ + ^ [n{pEp) -n{-pEp)] |. 
Finally, the comparison of the matrix elements of both sides of this equation yields, 

^P""^^^[^^^"('^^P)]' S^P = -^[l + 2«(;6£p)], (3.5) 



where we have used equations <3.H . ( 13. 2t and the evident property n[-x) = -1 - n{x). Here, 



^p = ep-Ai, £p=Jf2_|^p|2 (3,6) 



and £p, Ap are determined by the energy functional according to equations d3.3> . 
In the case of a homogeneous system, equation )2.18> is reduced to 

-^-M^o^O. (3.7) 

Equations B.SM3.7) provide a complete description of a homogeneous Bose system with single-particle 
and pair condensates if the energy functional £'(/'^,i/>) is given. The quantities bo and Ap should be 
interpreted as the order parameters associated with single-particle and pair condensates, respectively, 
and Ep as the quasiparticle energy. 

Let us now obtain the system of coupled equations using an explicit form of the energy functional 
given by equations <2.13> . <2.14) . In this case, equation I l3.7t for the condensate amplitudes takes the form 

bo [novm - m] + V I^/p [^(0^ + ^(P)] + 17 E gp^(P) = 0' (3.8) 



where we have used equations and introduced the single-particle condensate density no -b^bol V. 
Next, using equations d3.3> . we can find the quantities £p and Ap through /p - /p * and gp - g'L^. The 
subsequent substitution of equilibrium correlation functions <3.5> gives equations for £,p- e^- and Ap . 
We should also eliminate the correlation functions in equation <3.8> . After some algebraic manipulations, 
we come to the desired system of coupled equations for determining i,^, Ap, bo'- 

" - ^ + "0 [^(°^ + ^^P'] - ;T7 E [^^0) + V (p' - p)] 1 1 - [ 1 + 2n[pEp,)] 
Ap = v(p)fo2 _ _L ^ V (p + p') ^ [ 1 + 2n iPEp,)] , 

foo[2?zov(0)-<fo-^Ao] = 0. (3.9) 

Note that these equations, as well as equation <3.8> . are invariant with respect to the following transfor- 
mations: bo ^ b'f^- boe^f, ^ ^'p- ^p, Ap ^ Ap = ApC^^f. The same system of the coupled equations 
has been also derived within other approaches |14, 58]. It has three types of solutions. The first one, 
with «o = and Ap - 0, describes the state with no broken symmetry. In this normal state there is neither 
single-particle nor pair condensate. The second type of solution, with rio-O and Ap 0, breaks U(.l) sym- 
metry and corresponds to the state with a pair condensate, which is similar to a condensate of Cooper 
pairs in the theory of superconductivity. This condensation has been studied by a number of authors 1 2C|- 
I22II59I] . Finally, the third kind of solution, with noi^Q and Ap ^ 0, characterizes the state with broken U (1) 
sjmtmetry containing both single-particle and pair condensates. Note that the derived equations have no 
solution of the type Ap -Q,no + 0. 
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The chemical potential is related to the total particle number. From equations <2.15) . <3.5K we find 
that the total particle density can be written in the form 

n=no + «pair + Wquas , (3.10) 

where 

"pair = ^E(^-i)[1 + 2«(/^£p)], (3.11) 




(3.12) 



Here, we have introduced the particle density in a pair condensate (pair condensate density) npair and 
quasiparticle density «quas- In the normal phase, i.e., in the absence of both condensates, the particle 
number coincides with the quasiparticle number, as in the normal Fermi liquid. In a superfluid phase, 
the particle number is always greater than quasiparticle number. 



4. Zero temperature and contact interaction 

The coupled equations 13. 9t can be further simplified in the case of zero temperature and contact 
interaction, v(p) - v- const. Indeed, in this case, there are no quasiparticle excitations, n [pEp] - 0, and 
aU particles are in the single-particle or pair condensates. Moreover, due to contact interaction, Ap does 
not depend on momentum, Ap = A, and the structure of equations allows us to consider bo and A as real 
and positive quantities. After simple transformations, equations <3.9) are reduced to 

^E(|-i1-M.A = 0, (4.1) 



A 

where 



p 

1+ — y- 



- vno = 0, (4.2) 



p2 

f„ - + a, a- 2vno - A. (4.3) 
*^ 2m 



The spectrum of elementary excitations [see equations <3.6H takes the form 



As we see, it has a gap 6 = ijp=o given by 

5= Va'^-A^^ \/4vno(vno- A). (4.5) 

We will address the issue of single-particle excitation spectrum and its behavior in the region of small 
momenta herein below. 

Equations 14. It . I l4.2t allow one to find hq and A as functions of chemical potential [or taking into 
account equations ( I3.10> . 1 13. 11> as functions of total density]. As we have already mentioned, besides the 
solution corresponding to the normal state, these equations describe the state with a pair condensate 
only as well as the state with both condensates. Equation <4.2) shows that the solution corresponding 
to the state with a pair condensate (no - 0, A 0) exists only in the case of attractive interaction, v < 
0. However, this state was found to be thermodynamically unstable fsill . Therefore, we will study the 
solution of equations 1 14. It . 1 14. 2> that describes the state with both condensates (no ?i 0, A ?t 0). Moreover, 
the interaction is assumed to be repulsive, v > 0. 

The total particle density at T = 0, according to equations <3.10) . <3.11K becomes 

n= no + «pair= «o + T77l^[7r-l]- (4-6) 

2V p \tp J 
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The comparison of this formula with equation iA.l) yields the expression for the chemical potential, 

Id^ A-2v{no - n). (4.7) 



(4.8) 



From equation i3.5\ we can also find the normal and anomalous correlation functions at T = 0: 



Jn 



2En 



A 
'2K 



The ground state energy is obtained by setting T = in equations <2.13t . <3.5) . In addition, we should 
take into account the conditions for spatial homogeneity US. It and equations <4.6K ( 14. 8t . The ground state 
energy density S^^^ - / V is found to be 



»(0) 



1 

■ v[n-nof + nA-—+—Y^{Ep-ip). 



(4.9) 



When deriving this result we have eliminated the terms Xp (l^^^p) ^^d p^/2m using equations <4.2> . 1 14. 3> . 
respectively. 

Let us obtain an expUcit expression for pressure P - -Q./ pV. Here, Q is the equilibrium thermody- 
namic potential which is found by substituting the equilibrium values of the correlation functions and 
condensate amplitudes into equation <2.16t . Using the fact that S = at T = 0, one obtains P - 
and, consequently, 

a2 1 



(4.10) 



Here, as well as in equation I l4.9t . the quantity is given by equation <4.3) . From equation <4.10> we can 
conclude that the pressure is positive when v > that is one of the stability conditions. 

We now return to the coupled equations <4.1> . 1 14. 2> . First, let us eliminate the chemical potential in 
equation 14. 1> using equation <4.7> . Then, we replace the summation by integration over the variable 
X - p^/2m in both equations. The integral that appears in equation i4.2\ diverges at the upper limit. The 
physical reason for this divergence is that we have taken a delta-like contact potential with zero radius 
of interparticle interaction. In order to remove the divergence, we introduce a typical length scale kq for 
the range of interaction. In subsequent numerical computations, we will assume that ro coincides with 
the value of a repulsive core of the interaction potential, which is tjrpicaUy equal to a few angstroms. 
When integrating over the variable x having the dimension of energy, we will cut off the divergence 
by xq - {hko)^/2m - h^/2mrQ. In addition, for numerical analysis of these equations it is convenient to 
introduce the following dimensionless quantities: x = x/xq, A = A/xq, ct = u/xq, fio = «o''o' ^ - "''0 ^^'^ 
consider a - AgfiQ - A [see equations <4.3H as the sought quantity, instead of «o- Consequently, equations 
<4.1> , 14. 2> in dimensionless form are written as 



a-i- A 



4^ 



OC 

— f 



dx 



x+a 



\/(x-i-a)2-A2 



1 + 



2;r2 J 



dx- 



where 



g-- 



\/(x-i-a)2-A2 
vm \/2xom vm 



1, 



h^ro 



(4.11) 



(4.12) 



(4.13) 



is a dimensionless coupling constant and v is related to s-wave scattering length a by the well-known 
expression v - Anh^alm. We see that the dimensionless coupling constant is determined by the ratio of 
the scattering length to the radius of the interaction potential. 



a 

g^Aji—. 
ro 



(4.14) 
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Equation <4.1H reflects the fact that the total particle density is the sum of the particle densities in the 
single-particle and pair condensates (the first and second terms, respectively). The coupled equations 
1 14. lit . 1 14. 12> are characterized by two free dimensionless parameters — the total particle density n and 
coupling constant g. Therefore, having specified these quantities, we can find a and A and thereby deter- 
mine the fractions of single-particle and pair condensates as well as other characteristics of the system. 

We now present the dimensionless expressions for the quasiparticle energy, pressure, and density of 
the ground state energy. According to equation d4.4> . the single-particle excitation spectrum in terms of 
dimensionless quantities has the form 

Ep = ^J{p + af-A2^ (4.15) 

where Ep - Ep/xo and p - pro/H. In order to find an appropriate expression for pressure, let us replace 
the summation by integration over the variable x - p^/2m in equation <4.10t . The divergent integral, 
as above, is cut off by xq. Moreover, along with the introduced dimensionless quantities, we define the 
dimensionless value of pressure P - PIPq where Pq - h^/mrQ. Hence, 

A^ 1 

P^g[n^-nl) + — + j^Jo{&,A), (4.16) 

where 

1 

Jo [d, a) = j dx 



In a similar manner, one obtains the ground state energy density in a dimensionless form, 

S^°^^]-nA-^+g[n-hof- -^Jo [d, A) , (4.17) 
2 8g l&n'^ 

where =^"»/Po- 

In conclusion of this section, we estimate the value of Pq for a superfluid ''He. The interaction poten- 
tial of this system has a strong short-range repulsion whose radius (core) is tq - 2.55 A. The atomic mass 
of helium is m = 6.65- 10"^"* g. For these values, we find Pq =; 15.3 atm which is somewhat less than the 
value of crystalUzation pressure ~ 25 atm. 



x+d - \/{x+d)^-A'^ 



5. Weak interaction 

Here, we study the asymptotic solution of equations <4.11> . d4.12t in the case of a weak interaction. 
These equations contain two parameters g/n^ and hg, whose values are assumed to be small. 



71 



g 

— « 1, ng« 1. 



As one can see from equations 1 14. 11> . <4.12> . d-A-0 when g-0 and, consequently, a =; A « 1 at small 
g. Therefore, from equation 1 14. 12> , one obtains 

A=a;[l-^[v^l-H2fi- \/2a]j, 

with an accuracy of the terms of the order of d^'^: 

A=i5t[l--^(l- V^]]. (5.1) 

Next, since d and A are small, we neglect the terms d^ and A^ in equation <4.1H . Then, after its integration 
and subsequent substitution of equation <5.1> . we come to the following equation for a: 

4\/2 .3/2 
371^ 



Ahg- 2d ^d+ — Trga 
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Its approximate solution at small a is 



a-2ng+ — T 



- -,3/2 



ng- - [ng) 



After the substitution of equation <5.2) into <5.1> we can obtain A with the same level of accuracy: 



3 [ng] 



3/2 



(5.2) 



(5.3) 



The particle density in a single-particle condensate is expressed through a and A as follows: ho - 
[a + a) /4g [see equation <4.11H . Therefore, the asymptotic solution given by equations i5.2\ . <5.3) allows 
us to find the particle number densities in both condensates: 



1 (-j^5l2 
3^2 V-*; . "pair - 3^2 



no^n- — {ng) , 



i I- ^3/2 



(5.4) 



As we see, the terms ~ {fig)^^^ account for the presence of a pair condensate. 

The next step is to calculate the ground state energy density and pressure in the case of a weak inter- 
action. Both quantities, according to equation d4.16t . <4.17K are determined by the integral Jo [d, A) which 
at small a and A has the following asymptotic behavior: 

r ~\ 9 16 \/2 r-i-, 

Jo[a,A)~Jo[a) =j - 

15 



(5.5) 



or taking into account equation <5.2K 



Jo [a, A) Jo [a] -^4{ngf - ^ [ngf^ 

15 



Therefore, using equations <5.2M5.4> . we come to the following dimensionless expression for pressure: 



n g 1 , ^2 4 , .K/a 



It is also easy to find the pressure in dimensional form. 



1- 



1 

27r2 



vm 16 [vm] ' n 



3/2„l/2 



(5.6) 



(5.7) 



^h^ro 5 

Having obtained an explicit dependence of pressure on density, we can write down the speed of sound u: 

-) 1 dP vn 



m dn m 



1- 



(5.8) 



The second term in equations iS.H . <5.8> accounts for the finite range of interaction potential, while the 
third term is a correction responsible for the presence of a pair condensate. The ground state energy 
density <4.17> is determined by the same integral Jo [a, A). Therefore, taking into account equations d5.3> - 
d5.5> we have 



5(0) 1^ g ^ f- ^2 8 r- -,5/2 

' 7 ng) + 7 ng) 



or in a dimensional form. 



AO) . 



vn 



1- 



vm 16 (mv)^'^n^'^ 
■ -I- ■ 



The ground state energy - <f V can be expressed in terms of the scattering length a, 

27iah^N^ 



mV 



2 a 128 la^N'^'^ 



1 + ■ 



71 ro 15 s/n \ V 



(5.9) 



(5.10) 



(5.11) 
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This formula almost coincides with those obtained by Lee and Yang fso]. The difference is in the second 
term which takes into account the final range of interaction potential. However, the renormalization of 
scattering length a[l- lalnra) in equation <5.11> gives their original result. 

Finally, we note that a ^ A in the limit v ^ and, consequently, a gap 5 in the quasiparticle spectrum 
<4.4) tends to zero. In this case we come to the Bogolyubov spectrum l^. 




In the same limit, the density of a single-particle condensate tends to the total particle density (the pair 
condensate density tends to zero). 



6. Small density and strong Interaction 

Here, we study the system of small density without any restrictions on the value of interparticle in- 
teraction. To this end, it is convenient to introduce new dimensionless variables x-Kia and 77 - alAfig 
instead of a and A. In virtue of their definition, these variables meet the inequalities < ^,77 < 1. Then, 
equations <4.11) . <4.12) in terms of x and 77 read 

;?(l + x) + 4("^)'"'?'"/i(l) = l, (6.1) 
--^h{4ngnAngnx) = h (6.2) 



where 



00 

/i (l) = j dx\/x 





1 



x+ 1 



(6.3) 



Jzl^ngqAngnx] = fdx ^ (6.4) 

J{x + 4ngr]f-{4ngnxf 



We remind the reader that the left-hand side of equation <6.1) is just the sum of the single-particle con- 
densate fraction no In and pair condensate fraction Wpair/i. respectively. 

Now, assuming that the density is small we do not make any restrictions on the value of interaction, 

ng« 1, 4 ("^)'" « 1- 

Then, from <6.H - t6r4t . we find the equations for determining x and rj in the zeroth order in small param- 
eters: 

77(1 + X)-1 = 0, l--|_i.o, 

whence 

1 _ 1 + [gln^] 

The explicit expressions for x arid rj give a - Aghrj and A - Aghxri. Note that in this zero-order approx- 
imation, the total particle density coincides with the particle density in the single-particle condensate, 
n- no (there is no pair condensate). 

Let us now obtain the explicit form for the ground state energy and pressure. Both quantities, accord- 
ing to equations <4.16> . 1 14. 17> . are determined by the integral Jo (a, A), which at small density is written 
as follows: 

Jo [a, A) = Jo [AfigT], AngT)x) ~ (Ahgqxf, 
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where % and ^ are given by equations US. St . It is easy to find that in the given approximation, the ground 
state energy density coincides with the pressure, 



or in a dimensional form. 



2ro{nh)^ + vm 

The obtained expression for pressure allows us to calculate the speed of sound: 



n 1 dP vnro(7ih) 
u - — - 



2 



mdn mrQ(TihY + vm^ 12 

We see that >Q for repulsive interaction. This fact indicates the thermodynamic stability of the system. 
Next, from equation <4.5) we find the dimensionless expression for the energy gap in the single-particle 
excitation spectrum. 



5= \/a^-h^^Ang\ — -. (6.6) 

V 2 + g/n^ 

In virtue of the small value of hg, this gap remains small even at quite strong interaction. 

In conclusion of this section, let us calculate the pair condensate fraction [the second term in equation 
<6.1H . In the leading non-vanishing approximation over the small parameter, it has the form 



'pair 



AngY"v'^'Ji[x]- 



Taking into account that Ji [x] =; [tiIA)x^, where j is given by the first formula in equations d6.5> . the pair 
condensate fraction is found to be 

"pair _ .1/2 [glTt'^f^ 



- — n 



4 \l + [g/n^)]"'[2 + [g/n^)f' 



It is easy to show that the obtained function increases in the range of small values of g and decreases for 
"strong" interaction. Its maximum corresponds to the point g^ =; 50.9 and does not depend on the total 
particle density in the approximation under consideration. This fact agrees with numerical results shown 
in figure [T] 



7. Numerical results 

7.1. Dilute gases 

Bose-Einstein condensation in dilute ultracold gases of alkali-metal atoms was first realized experi- 
mentally in 1995 |l-3]. The particlenumber density in the condensed atomic cloudis n ~ 10^^ + 10^^ cm~^. 
Then, the corresponding dimensionless density in equations <6.1> . <6.2> is a small quantity, n ~ 10"^^ + 
10~^, where we have taken ro = 3 A. The binary atomic interaction in such systems is usually approx- 
imated by the contact interaction potential expressed through the scattering length a. For example, for 
^^Rb and ^^Na, the scattering lengths are equal to a =; 90ao and a =; 19.1flo. respectively, where ao =; 0.53 A 
is the Bohr radius 1 6]. From equation <4.14t , we can find that the dimensionless coupling constant: g^ 200 
for ^''Rb and 42 for ^^Na. The numerical analysis of equations I IB.H . I |6.2> for a system with parameters 
of dilute gases of alkali-metal atoms are presented in figures[l][2] 

Figure [1] shows the dependencies of the pair condensate fraction fipairin on the coupling constant g 
for the above shown values of the total particle density n (all quantities are dimensionless). The curves 
have nonmonotonous character. Moreover, their maximum is reached at g^ =; 50 and practically does 
not depend on the total density. This result, as we have already seen in the previous section, follows 



13603-13 



A.S. Peletminskii, S.V. Peletminskii, Yu.M. Poluektov 




from the analytical study of equations dS.lt . <6.2) . We can conclude that the pair condensate fraction is 
several orders less than the single-particle condensate fraction and, consequently, the dilute systems are 
described by the Gross-Pitaevskii equation with a high level of accuracy. In addition, in the limit 0, 
the pair condensate density also tends to zero and thereby no h. 

Figure [2] presents the dimensionless energy gap 8 in the spectrum of single-particle excitations as a 
function of dimensionless coupling constant g for various values of the total particle density h. At small 
densities, the gap remains small even for large values of the coupling constant. The gap tends to zero in 
the Umit ► which agrees with Bogolyubov's result 

7.2. Model of superfluid ^He 

Superfluid ^He is a strongly interacting system of sufficiently high density in which Bose-Einstein 
condensation is also manifested |48-S1]. The density of a liquid helium at zero temperature and low 
pressure is about n ~ 2.18- 10^^ cm~^. The atomic mass is m = 6.65- 10"^^ g. The interaction potential of 
helium has a strong short-range repulsion whose radius is ro ~ 2.55 A, followed by a weak intermediate 
range attraction. Therefore, the dimensionless density n - nr^ in equations 16. 1> . <6.2> is n =; 0.36. The 
scattering length a for *He atoms varies from 46.1 A to 100 A depending on the interaction potential |l6l|l. 
Hence, equation d4.14> shows that the dimensionless coupling constant g lies in the interval from 227 to 
492. The numerical analysis of equations <6.1) . <6.2) for a superfluid ^He are presented in figures[3l|4l 




Figure 3. Pair condensate fraction Kpair/w as a func- 
tion of the coupling constant ^at the density of liquid 
*He, h = 0.36. 




Figure 4. Energy gap 5 in the single-particle excita- 
tion spectrum as a function of the coupling constant 
gat the density of liquid *He, n = 0.36. 
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Figure[3]illustrates the behavior of the pair condensate fraction depending on the coupling constant. 
We can see that it significantly exceeds the single-particle condensate fraction. For the above shown val- 
ues of the coupling constant, the single-particle condensate fraction is just about 5% + 10%. This agrees 
with the condensate fraction in a superfluid ^He measured in experiments fS-^Sl*] and found within 
Monte Carlo calculations [S21 and theoretical predictions based on the computation of single-particle 
density matrix which is expressed through experimentally measured structure factor (s^lsjl (see also 
reference fss'l and references therein). 

FigurelDpresents the dependence of the dimensionless energy gap S in the single-particle excitation 
spectrum on g at helium density. For example, at g- 227 we have 8 =; 13.45 or in temperature units 
d =: 12.5 K. This value is somewhat greater than the value of the roton gap Arot ~ 8.65 K and less than the 
maxon energy Amax ~ 14 K. The Landau gapless spectrum can be a result of superposition of excitations 
of various nature. Moreover, some experiments (g^ directed toward studying the excitations in a 
superfluid ^He also show the complex structure of Landau spectrum. 

8. Discussion 

Within the quasiparticle approach based on the maximum entropy principle, we have given a deriva- 
tion of the coupled equations describing a superfluid Bose system with single-particle and pair conden- 
sates. These equations have been analyzed both analytically and numerically. In order to remove the di- 
vergences, we introduce a typical length scale ro for the range of interaction potential. For dUute systems 
(alkali atomic gases), the pair condensate fraction is negligibly small at zero temperature, «pair/i ~ 10"^ 
and it grows with an increasing density n. For a superfluid *He, the role of pair condensate at r = is 
dominant because the single-particle condensate fraction no/nis less than 10%, which agrees with experi- 
mental data on neutron scattering in a superfluid heUu m kSHSlll . with Monte Carlo calculations 1 52], and 
with other theoretical predictions (see references jssMsal and references therein). We have also found 
other characteristics of the system such as pressure, ground state energy, speed of sound, and single- 
particle excitation spectrum. Note that the latter exhibits a gap. We believe that this fact is sufficiently 
justified because the used approximation introduces quasiparticles by the general way and respects the 
basic principles of statistical physics. The studied approach can also be used to derive equations of two- 
fluid hydrodynamics predicting the second sound wave [41, 46]. As we have already mentioned in the 
introduction, the analysis of Bogolyubov's -theorem does not provide a general conclusion concern- 
ing the excitation spectrum of a superfluid, and long-wavelength density excitations are insensitive to the 
gauge symmetry breakirtg [28] . The latter fact is directly confirmed by the neutron scattering experiments 
in a superfluid *He je^lil]. The single-particle excitation spectrum exhibits a gap in the presence of a 
pair condensate because the pairs have the dissociation energy and, consequently, in addition to phonon 
branch of spectrum, there exists another branch corresponding to the excitation of pairs [29,1. The same 
qualitative structure of the spectrum is obtained within the proposed approach. Note that the existence 
of a gap resembles the situation in a solid: if there is more than one atom per unit cell, both acoustic 
and optical branches appear. Moreover, as we know, a neutral Fermi superfluid is characterized by the 
single-particle excitations with a gap as well as gapless phonon mode. The separation of single-particle 
and collective excitations is probably less marked in a Bose system than in a Fermi system as a conse- 
quence of "hybridization" of these branches due to the presence of a single-particle condensate |64, 6Sj]. 
Therefore, one can expect that the Landau spectrum is a result of superposition of excitations of various 
nature. It is worth stressing that the experiments (gl, directed toward studying the excitations in a 
superfluid *He also show its complex structure. 
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Po/ib OAHonacTi/iHKOBoro ra napHoro KOHAeHcaTie y 
6o3e-ci/iCTeMax 13 AoembHOio iHTeHCi/iBHicTK) B3aeM0Ai'i 

O.C. ne.neTMi/iHCbKi/ibi, C.B. ne.neTMi/iHCbKi/ibi, fO.M. Flo/iyeKTOB 

iHCTHTyr TeopeTH4HoT (fiiswKH Im. O.I. Axiesepa, HHI4 XttTI, ey/i. AKafleMi^Ha, 1, 61 108 XapKie, YKpaiHa 

BuBHaeibCfl Hafln/inHHa 6o3e-cwcTeMa 3 oflHO^acTUHKOBUM ra napHUM KOHfleHcaraMw Ha ochobI naniBcfieHO- 
MeHoyiori4HoT reopii 6o3e-piflWHW, UKa He npwnycKae c;ia6K0CTi MixHacTWHKOBOi B3aeM0flii. I3 BapiaL4itiHoro 
npnHL4nny fl/i^i eHiponiT BWBefleno cucreMy piBH^Hb, 1140 onwcyKJib piBHOBaxHwti crah laKoT cucreMH. U,\ piB- 
HflHH^i npoaHa;ii30BaHO y BwnaflKy hy/ibOBoT leMnepaiypn aHa/iiTW4H0, laK i Huce/ibho. noKa3aHO, mo 4ac- 
TKH OflHOHacTUHKOBOro Ta napHoro KOHflencaTiB cyrreBO 3a;iexaTb Bifl noBHOi rycTHHH cucreMH. flpn rycTunax, 
mo flocaraiOTbCfl y KOHfleHcarax aroMlB /lyxHwx Mera/iiB, MaPixe Bci HacTWHKW 3HaxoflflTbCfl b OflHO^acTUHKO- 
BOMy KOHfleHcaTi. flOBefleHO, mo npn rycTHHi piflKoro re/ilK) nacTKa OflHonacTUHKOBOro KOHflencaTy e MeHLUOK) 
3a 10%, mo ysroflxycTbca 3 eKcnepuMeHia/ibHi/iMn flaHUMM 3 HenpyxHoro HetiTpoHHoro po3ciK)BaHHfl, po3pa- 
xyHKaMM MBTOflOM MoHie-Kap/io Ta Ihujumw TeopernHHUMn nepefl6aHeHHflMW. BHaPifleHO BMpa3M fl/ia eHeprii 
ocHOBHOro craHy, rwcKy, cthc/iuboctI cwcreMU, mo BUBHaerbcn. npoaHa;ii30BaHO laKOx cneKip OflHonacTHH- 
KOBI/IX 36yflxeHb. 

K;iioMOBi c;ioBa: Haflnjii/iHHicTb, 6o3e-eiiHinTeiiHiBCbKa KOHfleHcai^ip, oflHOHacji/iHKOBi/iii ra napnuii 
KOH^eHcam, KBaai^acTi/iHKi/i, cneKrp 36yA>KeHb 
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